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Abstract. We present the complete next-to-next-to-leading-order calculation of the correlation function of two Polyakov 
loops for temperatures smaller than the inverse distance between the loops and larger than the Coulomb potential. We discuss 
the relationship of this correlator with the singlet and octet potentials which we obtain in an Effective Field Theory framework 
based on finite-temperature potential Non-Relativistic QCD, showing that the Polyakov loop correlator can be re-expressed, 
at the leading order in a multipole expansion, as a sum of singlet and octet contributions. We also revisit the calculation of the 
expectation value of the Polyakov loop at next-to-next-to-leading order. 
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INTRODUCTION 

The Polyakov loop and the correlator of two Polyakov 
loops are the order parameters of the deconfinement tran- 
sition in pure gauge SU(N) theory As such, they 
have been extensively measured on the lattice, both in 
pure gauge and in full QCD (for a review see [3]). From 
a phenomenological point of view, the Polyakov loop and 
the correlator are important because they are related re- 
spectively to the free energy of a static quark and of a 
static quark-antiquark pair in a thermal bath. As such, 
they can be used as input for the phenomenology of 
heavy quarkonia in the quark-gluon plasma through po- 
tential models (see for instance 140). However, although 
the spectral decomposition of the correlator is known 
its relationship to the in-medium heavy quark po- 
tential is still in need of clarifications [7], especially after 
the recent developments in a derivation of the potential 
from QCD [8, 9] and its rigorous definition and calcula- 
tion within an Effective Field Theory (EFT) framework 

The work of Ref. 11311 . which we shortly present here, 
aims at a better understanding of the Polyakov loop cor- 
relator in the short distance regime rT <C 1. This is 
achieved by performing first a perturbative NNLO (order 
g 6 ) computation of the correlator, improving the leading- 
order results of 10, [Till and complementing the order g 6 
results of Nadkarni IT15ll , which are valid for rmo ~ 1 , mo 
being the Debye screening mass. As a side result, we also 
obtain the perturbative expression of the Polyakov loop 
at the NNLO (order g A ). We then proceed to construct 
an EFT framework, based on potential Non-Relativistic 
QCD (pNRQCD), which allows us to reobtain the pre- 
vious result for the correlator in terms of the colour- 



singlet and colour-octet zero-temperature potentials of 
pNRQCD and further contributions, which we compute, 
induced by the thermodynamic al scales. 

PERTURBATIVE RESULTS FOR THE 
POLYAKOV LOOP AND ITS 
CORRELATOR 

We briefly present here the calculation methodology 
and the results for the Polyakov loop and its corre- 
lator. The former is defined as (Lr) — (TtLr), where 

Lr = Pexp Mg f^ T dtA°(x, t)J is the Polyakov line in 

a colour representation R, either the fundamental F or 
the adjoint A, and Tr is the colour trace normalized by 
the dimension of the representation R. The correlator is 
defined as (frl* (0)frL F (r)); what we consider in the 
following is its connected part Cpl(^) in the fundamental 
representation. 

We adopt the static gauge lfl7tl . defined by the condition 
d°A° = 0. In this gauge the Polyakov line becomes lo- 
cal in time, i.e. L = exp ((igA°(x))/T) . The one-loop ex- 
pression of the 00 component of the gluon self-energy 
in this gauge, which is a fundamental ingredient in both 
computations, was not fully known in the literature and 
has been computed in 11311 . 

Both calculations assume a weak-coupling hierarchy 
T ^> mo, where m D — g 2 T 2 /3(N + nj/2), with rif light 
quarks. In the case of the correlator we furthermore as- 
sume 1 jr ^s> T ^> mo g 2 1 r, where the last scale is the 
Coulomb scale that appears naturally in a perturbative 
quark-antiquark correlator. We stop our perturbative cal- 
culation at order g 6 (rT)°, meaning that we do not show 



terms with powers of g greater than 6, and at order g 6 
we do not show powers of rT greater than zero. The 
aforementioned hierarchies are implemented by separat- 
ing the momentum regions in the integrations by appro- 
priate expansions, and by resumming the Debye mass in 
the longitudinal propagators when the loop momentum 
is of the same order. All divergences are treated in di- 
mensional regularization; the Polyakov loop is finite at 
the order considered, whereas the correlator requires the 
usual charge renormalization of QCD. 
Our result for the Polyakov loop is 
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Our result differs from the long-time accepted result of 
Gava and Jengo [16] in the finite terms at order g 4 and 
agrees with another recent determination lfl8ll . For the 
correlator we have 
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where the scale of the running coupling in the first, 
leading-order term is fixed to 1 jr by the higher-order 
terms. Nadkarni [ 15] computed the same quantity under 
the assumption T ^> l/r ~ nip; as such, a direct com- 
parison between the two computations is not possible. 
However the part of our results that does not depend on 
the hierarchy rT <C 1 agrees with Nadkarni's results ex- 
panded for rniD ^ 1 . 



THE POLYAKOV-LOOP CORRELATOR 
IN AN EFT LANGUAGE 

We now proceed to briefly illustrate the principles and 
the construction of the EFT framework of foil . Such a 
framework allows us to reobtain the result of Eq. (f2]l in 
terms of colour singlet and colour octet correlators, mak- 
ing more explicit the physical meaning of the calculation. 
The basic principle is that we implement the hierarchy 
l/r ^> T ^> mo ^> g 2 /r at the Lagrangian level, by in- 
tegrating out from (static) QCD these scales. We start 



from the largest, the inverse distance l/r. This leads to 
pNRQCD. Since 1 jr is larger than all other scales, ther- 
mal ones included, those can be put to zero in the match- 
ing and the Lagrangian is the usual zero-temperature one 
11911 . In Euclidean space-time the action reads 
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where f y = f d 3 y, the trace is over colour indices, S = 

S/VN and O = \20 a T a are the colour-singlet and 
colour-octet quark-antiquark fields and E' — F'° is the 
chromoelectric field. The fields S and O depend on the 
relative distance r, on the center-of-mass coordinate x 
and on the Euclidean time T. V A , Vb and Vc are the match- 
ing coefficients of the first operators in the multipole ex- 
pansion; for our purposes they are equal to unity. V s and 
V are the singlet and octet potentials in pNRQCD. The 
former is known to three loops l2(ill . whereas the latter 
to two loops 02111 . For our purposes the two-loop result 
l22tl is sufficient. Finally 5Jz? p nrqcd contains operators 
in the multipole expansion of order higher or equal to r 3 , 
which are shown not to contribute to the Polyakov loop 
correlator at the intended accuracy. 
We are then able to write the connected Polyakov loop 
correlator in terms of singlet and octet correlators as 
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The suppressed terms in the last line originate from the 
8 Jzf pnrqcd terms in the Lagrangian. We have thus shown 
that at the leading order in the multipole expansion of 
pNRQCD the Polyakov loop correlator can be expressed 
as a sum of singlet and octet correlators. The matching 
at the scale l/r furthermore yields Z s =7^ = 1, consis- 
tently with the spectral decomposition of the Polyakov 
loop correlator j5l|6J], whereas the (SS^) and (00^) cor- 
relators read 

(SS^) = e- V *W T {\ + 4) = e S^ T ^l\ (5) 
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where we have defined the gauge-invariant singlet and 
octet free energies f s and f from these correlators, (La) 
is the Polyakov loop in the adjoint representation and 8 S 
and 8 () stand for the loop corrections to these correlators. 
(Lf) (La), 8 s and 8 encode thus the contribution of the 
lower scales, in particular T and mo, that are still dynam- 
ical in the EFT. We can proceed to compute 8 S and 8 by 
computing loop diagrams in pNRQCD, again separating 
the contribution of T from the one of mo with the previ- 
ously exposed techniques. In the end, plugging back the 
results for 8 S , 8 (see II 1310 and those for the Polyakov 
loops (see Eq. (Q])) in Eq. ©, large cancellation happen 
between those terms and we reobtain the result for Cpl (r) 
previously shown in Eq. (|2), which was obtained from a 
direct perturbative computation. 

In comparison with the recent developments in the 
construction of a real-time finite-temperature pNRQCD 
framework for the study of in-medium heavy quarko- 
nium [fToL [l2tl . we observe that the free energies f s and 
f do not agree completely with the real parts of the real- 
time potentials therein computed. The difference can be 
traced back to the different boundary conditions: cyclic 
imaginary time in this case, large real time in that case. 
EFT-based calculations of the correlator were performed 
in the past for different scale hierarchies; beside the al- 
ready mentioned work of Nadkarni for \/r ~ mo in 
EQCD [15], we also have a calculation within MQCD 
in the magnetic screening region mo 3> 1/r II23I1 . In both 
cases the scale \/r was however not integrated out and 
the complexity of the bound-state dynamics remained 
implicit in the correlator. 

Finally in Ref. Q5|] the gauge structure and transformation 
properties of allowed intermediate states contributing to 
the correlator of Polyakov loops were analyzed; we ob- 
serve that our results are in agreement with this analysis, 
since in the weak-coupling regime we considered, both 
our colour singlet and colour octets transform accord- 
ing to the allowed transformation; for strongly-coupled 
pNRQCD the allowed intermediate states would be the 
colour singlet and its hybrid gluonic excitations. 
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